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Abstract 

The distribution of collisionless particles with infinite motion in the presence of 
a local potential well is discussed. Such distribution is important for interpretation 
of results of dark matter searches. The relationship n/v=const, where n and v are 
respectively number density and velocity of particles, is derived for particles crossing a 
local potential well. The limits of application of this relationship are specified. 



The question on the distribution of collisionless particles inside the local potential well 
is important for the problem of Dark Matter (DM) searches (e.g., ^P-JHl)- In the papers [5] 
a relationship for number density of DM particles n inside local potential well 

n ,„ . 

— = const, (1) 

v 

where v is the local velocity of DM particles in the reference frame related with the well, 
was used putting aside derivation. We adduce derivation of Eq. (JIJ in the present note and 
specify the limits of its application to the problems of dark matter search. 

In general, the particles which experience an infinite motion are of interest. They are 
supposed to be distributed homogeneously at sufficiently large distance from the potential 
well. Potential well is supposed to be spherically symmetric and small in size, so the particles 
cross it quickly and the effects of scattering can be neglected for the bulk of particles. These 
conditions, as a rule, arc fulfilled for WIMPs crossing the Solar system. 

Let us introduce necessary notations. Intensity of the flux 

j — dN (2) 

ds„dtdn' 1 ; 
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where N is the number of particles, S n is the area normal to given direction, £1 is the solid 
angle. If particles are distributed in velocity then dv should be added in denominator of 
Eq. J2J . In case of isotropic monochromatic distribution 

The number of particles passing through given area dS per unit time is 

dN = J IdS cos Odfl, (4) 

where 9 is the angle between direction of incident particles and normal of area. In case of 
isotropic flux, N from one side of area dS (0 < 8 < tt/2) is 

711} 

dN = irI-dS= —dS. (5) 

4 v ; 

For sphere of radius R it gives the well known expression 

nv 2 

N on sphere = ~ 4vri? = nV<7, (6) 

where a = ttR 2 is the cross section of sphere. 

Velocity of collisionless particles in a given point of potential well (R) can be defined 
through the relevant escape velocity v esc {R) 



vr 



= Vv 2 oo+V 2 ese(R)- (?) 



Throughout in the paper indexes "R" and "oo" correspond, respectively, to the points at 
the distance R and at "infinity" (at sufficiently large distance, at which the influence of 
potential well is negligible). 

Let us first consider a potential well of gravitating body with radius R (e.g., a star). For 
particles moving from infinity the cross section of their capture by this body is given by 



V w 5o / "So 

Here Co = ttR 2 is the geometrical cross section of the body, v esc = v esc (R) is the escape 
velocity on its surface, and the Eq.0 was used to derive the last expression in the Eq.JSJ). 
The number of particles, being isotropically distributed on infinity (though the following 
formula is not restricted by this case only), falling on star surface per unit time will be 

N = n 00 v 00 a oc . (9) 

From the other hand the same number can be expressed through the corresponding values 
near the surface. If we suppose that distribution near the surface is also isotropic, we obtain 

N = n R VR<T . (10) 

Equating Eas. (|9ll0|l one gets 

— = — . (11 

VR Woo 

This derivation gives a hint for the Eq.Q), requiring to resolve the question of (an)isotropy 
inside potential well and to generalize the expression on the case of more realistic distribu- 
tion. 

Let us consider particles on infinity, having arbitrary angular distribution of their velocity 
vector (voo). Assume that particles move in a central-symmetric potential field. The shape 
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Figure 1 : Illustration of motion of test particle in a central-symmetric attractive field with 
the center O. 



of potential well is not specified. It will be enough in this case to use angular momentum 
conservation. 

Let particles on infinity have the only value Voo = \voo\- One separates from them the 
flux concentrated along a given direction loodtt. Intensity can be factorized as 

loo = ™ooKoo/oo(#, <f>), (12) 

where /oo($, <P) is distribution function in angle coordinates of Woo relatively to the centre of 
attraction normalized on unit (e.g., for isotropic flux Eq.10 gives f x = l/(47r)). 

Impact parameter p, angle \& of incidence on sphere of radius R around gravitating center 
(equipotential layer) are attributed to the motion of test particles (with mass m) of given 
flux as shown on Figure 1. 

Angular momentum conservation law gives us 

tnVoop = mvnRsm^. (13) 

Position where the particle crosses the sphere C\ is not specified. Note that by force of 
Eq. l|13|l an incident angle is simply defined by R at the given Voap. In particular, an equality 
of the incident angles in the points C\ and C2 follows. For the further use one can express 
dp 2 through cos iff (not caring about the sign to be eliminated by proper choice of integration 
limits) 

\ 2 



dp =2 — cos I'd cos*. (14) 

The number of particles falling from the outside with the fixed impact parameter p...p + dp 
or corresponding incident angle * per unit time on sphere of radius R is given by 

dNoutside = IoodVt ■ 2npdp. (15) 

For the number of particles crossing the given sphere both from the outside and from the 
inside due to the above mentioned symmetry we have dN sp h — dN outS id e + dNi ns id e — 
2dN outside- Transforming it with the use of Ea. p4|l one obtains 

dN sph = I^dfl ■ AttR 2 (^—^j cos^dcos*. (16) 

One can notice from the comparison of Eq. 1)16(1 and Eq. J3J that intensity of particles falling 
on sphere (or near it) integrated over asimuth angle of vr (where angular coordinates are 

defined relative to R) is dlf;^) = IoodQ ■ (j^^j — const. 
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Consider particles within spherical layer of radius R. Let us calculate the number of 
particles present at the same time inside the layer R...R + dR. Each particle moving through 
the layer with incident angle W will reside in it during time 

dR 

dr = -. (17 

Vr cos w 

Then the mean number of particles present in the layer will be 

dN layer = dN sph dr = I^dVl • — — ( — \ dcos (18) 

VR \VooJ 

Note that the distribution on cos \& of number of particles (of their arbitrary flux from 
infinity) inside the layer is constant, i.e. 

= const, (19) 

dcos* 

It is worth to emphasize that the angle \& is defined relative to the direction towards the 
center of attraction (O) R wherever the point C\ is situated on the sphere. It does not imply 
space homogeneity of distribution of the points C\ on sphere. 

After trivial integration of Eq. (|18|l over cos W in the limits from to 1 (cos * = 
corresponds to maximal p — ^f^) one writes with the use of Eci. (|12fl 

delayer "oo r / r, ,\ J( -. /r»n\ 

-vr ■ foo(9,(j))dn. (20) 



4TrR 2 dR v. 



OC 



The value 4^2^ is the number density of particles averaged over sphere - equipotential 
layer, ur. After integration over dfl - directions of incident particles from infinity, one 
obtains 

= ^L. (21) 

VR Voo 

In case of isotropic distribution foo at infinity one can generalize Eq. 1)19(1 as an averment 
of isotropy of incident flux locally in any point of (any) equipotential layer. Ea. 1(21(1 in this 
case takes the form of Eas. (lll(l .l(TJl (i.e. averaging over the layer is not necessary). 

To generalize given above arguments for a case of arbitrary distribution in velocity 
Voo, one introduces dv in denominator of Eq.(0l. One separates the flux Ioodvoodfl = 
"■00^00/00(^00, 0, (p)dvoodn. Then En. 1(20(1 takes the form 

Ti 

dnR = —vr ■ foo(voo,0, (f))dVoadQ,. (22) 

Woo 

For any point i, a notation dni = riifi(vi, 9, (f>)dvid£l = riifi(vi)dvi, where indexes for angles 
are omitted and / fi(vt, 9, 4>)dQ, = fi(vi), is adopted (n* and fi are overlined, if any). Veloc- 
ities Voo and vr are mutually related by Eq.0 and for their differentials the relationship 

VRdVR = Doorffoo (23) 

is valid. Integration of Eq. 1(22(1 over d£l with introduced notation yields 

dnR = dupo ^ 

VR Voo 

So, account of distribution in velocity leads to replacement rii — > drii(vi) in the Ea. 1(21(1 . 

Velocity distribution, averaged over sphere of radius R, is determined through that on 
infinity as 



Mvr) = ^L^-f^VooivR)). (25) 

n>R vio(VR) 
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Its normalization on unit is provided by factor Note, that /r(«.r) has a threshold 

vr > v esc (R), where it experiences a jump from zero. 

The relations for number density considered above are applied to the particles with 
infinite motion and do not take into account particles captured by potential well. The 
latter would contribute in range just below the threshold v esc (R) (and in perfect case would 
reproduce equilibrium form of velocity distribution). In case of WIMPs trapped by Solar 
system, still there is no definite answer on the effectiveness of respective mechanisms of the 
capture and, especially, of the escape of WIMPs (see, e.g., [7], [5], @]). 

Application of Eq. Q or Eq. (|24|l for estimation of accumulation of WIMPs by the Earth 
can cause errors. The Earth is inside Solar potential, what could be taken into account with 
"two-step" application of given equation. One can get first the proper WIMPs parameters 
inside the potential of the Sun at the distance from it corresponding to the Earth orbit and 
then pass to the Earth reference frame and get those parameters inside the Earth potential. 
However, the use of Eq. l|24[l at first step will give density distribution averaged over all 
equipotential layer around the Sun rather then that along the Earth's orbit only. 

Because of the averaged form of the obtained result Eq.QJ can not be applied to the 
analysis of annual modulation in direct WIMP searches |10j . 

One of the implications of given here equations can be connected with predictions of 
possible existence of anti-meteorits in our Galaxy Analysis of the effects of their 

annihilation on Sun can provide crucial test of this prediction. For estimation of frequency 
of such events Eq. can be used and to study observable effects (specific forms of flares) 
produced by them Ea. H19f) can be useful. 

Generalization of Eq. Q on case of Z?-dimensional space has the form 

n ,„„. 
n _„ = const. (26) 

v 

In case D = 1 it leads to the familiar result from theory of gas flow in tube with variable 
cross section. 
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